Fungible dynamics: there are only two types of entangling multiple-qubit interactions 
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What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum 
system? It has been shown that all two-body Hamiltonian evolutions can be simulated using any 
fixed two-body entangling n-qubit Hamiltonian and fast local unitaries. By entangling we mean 
that every qubit is coupled to every other qubit, if not directly, then indirectly via intermediate 
qubits. We extend this study to the case where interactions may involve more than two qubits 
at a time. We find necessary and sufficient conditions for an arbitrary n-qubit Hamiltonian to be 
dynamically universal, that is, able to simulate any other Hamiltonian acting on n qubits, possibly 
in an inefficient manner. We prove that an entangling Hamiltonian is dynamically universal if and 
only if it contains at least one coupling term involving an even number of interacting qubits. For 
odd entangling Hamiltonians, i.e., Hamiltonians with couplings that involve only an odd number 
of qubits, we prove that dynamic universality is possible on an encoded set of n — 1 logical qubits. 
We further prove that an odd entangling Hamiltonian can simulate any other odd Hamiltonian and 
classify the algebras that such Hamiltonians generate. Thus, our results show that up to local unitary 
operations, there are only two fundamentally different types of entangling Hamiltonian on n qubits. 
We also demonstrate that, provided the number of qubits directly coupled by the Hamiltonian is 
bounded above by a constant, our techniques can be made efficient. 

PACS numbers: 03.65.Bz, 03.67.-a 



I. INTRODUCTION 

One of the main goals of quantum information science 
is to characterize the physical resources that are universal 
for quantum computation and simulation. Recently, the 
role of entangling quantum dynamics has been studied in 
depth, and it has been shown that a fixed two-body en- 
tangling Hamiltonian evolution acting on a system, plus 
the ability to intersperse local unitary operations, can be 
used to efficiently simulate any other two-body Hamil- 
tonian and hence is universal for quantum computation 
(see, for example [10HSI!,lllEEI!,[l3, 

and refer- 
ences therein.) 

We extend this study to Hamiltonians containing in- 
teraction terms involving more than two qubits at a time. 
Using results from quantum control theory, we determine 
what dynamics can be simulated with an arbitrary fixed 
rt-qubit Hamiltonian, and complete local control in the 
form of one-qubit unitary operations. In particular, we 
derive a simple criterion determining which Hamiltonians 
are universal given local unitary operations. 
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In contrast to two-body entangling Hamiltonians, 
which are always universal given local operations 0,07 0, 
000000liil. we W111 see that n °t all many-body en- 
tangling Hamiltonians are universal, given arbitrary one- 
qubit operations as the only additional resource. This 
was previously noted for a specific example by Vidal and 
Cirac |9j- Thus, unlike two-body Hamiltonians, many- 
qubit entangling Hamiltonians are not all equivalent up 
to local operations. The situation is somewhat analogous 
to the study of multi-party entangled states, where it has 
been shown that there exist different types of entangle- 
ment, inequivalent up to local operations and classical 
communication 0,0,0 (LOCC). For example, it is 
now well known that the "GHZ" and "W" states are not 
equivalent up to LOCC. 

Consider the following illustrative example. Suppose 
we are given a Hamiltonian acting on three qubits, H = 
X<g)X®I + I®X®X, where X is the usual Pauli a x 
operator. Then, if we can perform arbitrary local unitary 
operations on each of the qubits, it has been shown in 
Ql 13 that it is possible to simulate any other 
Hamiltonian interaction on three qubits, such as H' = 
We call such a Hamiltonian universal. On 
the other hand, given H' and arbitrary local unitaries, 
it turns out that it is not possible to simulate H, and 
thus H' is not universal. A proof of this is given in 0] ; 
both these results will also follow from the general results 
given in this paper. 
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We address the problem of universality in full gener- 
ality by giving a necessary and sufficient condition for a 
Hamiltonian to be universal. Our condition is dependent 
only on simple properties of the Hamiltonian's decom- 
position into tensor products of Pauli operators. The 
proof of this condition is constructive, in the sense that 
it provides, in principle, an algorithm for using a univer- 
sal Hamiltonian to simulate any other interaction. How- 
ever, the techniques used in our construction are not es- 
pecially practical, especially in the presence of noise, and 
it remains to be seen if more practical constructions are 
possible. 

In addition to our criterion for universality, we also 
examine what can be done when a given Hamiltonian is 
not universal. In particular, we show that there always 
exists a simple encoding scheme to make these Hamilto- 
nians universal. 

Let us make a more precise statement of our results. 
Suppose we are given a Hamiltonian H acting on n 
qubits, which can be written uniquely in terms of its 
Pauli operator expansion 
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H = zZ h h-i» a h ® ■ ■ ■ ® a ir, , (!) 

jl,-,jn=0 

where the hj 1 ,_j n are real numbers and a\, a%, 03 are the 
Pauli sigma matrices X, Y, Z, respectively, with 00 = I 
the identity. We say that a subset S of the qubits is 
coupled by this Hamiltonian if there is a non-zero term 
in H coupling those specific qubits. 

When is H universal? An obvious condition is that the 
set of qubits coupled by H must be connected. That is, 
it should not be possible to partition the qubits into non- 
trivial sets, S and S, such that every term in the Pauli 
operator expansion couples either a subset of S or a sub- 
set of S. If this were the case then H could not be used to 
generate entanglement between the qubits in S and the 
qubits S, and thus would not be universal. We say that 
a Hamiltonian connecting all the qubits in this way is an 
entangling Hamiltonian. Note that this definition may 
be restated in the language of graph theory: if qubits 
correspond to vertices in a hypergraph, and couplings 
between qubits correspond to hyper-edges, then the con- 
dition that the Hamiltonian is entangling corresponds to 
the condition that the hypergraph is connected. 

With this background, our main results are easily 
stated. Our first result is that an entangling Hamilto- 
nian and local unitaries is universal if and only if the 
Pauli operator expansion for the Hamiltonian contains 
a term coupling an even number of qubits. This result 
provides a simple, easily checkable criterion to determine 
whether or not a Hamiltonian is universal. Returning to 
our previous examples H = X®X®I + I®X®X and 
H 1 = X ® X ® X, this criterion tells us that H is uni- 
versal when assisted by local unitaries, while H' is not, 
in agreement with the earlier claims. 

Our second result concerns what happens when the 
Pauli operator expansion contains only odd terms, and 



thus is not universal. We say that such a Hamiltonian 
is an odd entangling Hamiltonian. We will prove that 
an odd entangling Hamiltonian acting on n qubits is ca- 
pable of simulating any other odd Hamiltonian acting 
on those qubits. Thus, the odd entangling Hamiltonians 
are a fungible physical resource, since having any one is 
equivalent to having any other, up to local unitary oper- 
ations. Furthermore, we show that an odd Hamiltonian, 
together with local unitaries, generates either the simple 
Lie algebra so(2 n ) or sp(2"), depending on the number 
of qubits n that are connected by the Hamiltonians. 

Our third result also concerns odd Hamiltonians. We 
prove that, with appropriate encoding, an odd entangling 
Hamiltonian and local unitaries is universal on a set of 
logical qubits. Thus, there is not too great a loss in 
space efficiency when one attempts to use such a Hamil- 
tonian to simulate an arbitrary interaction. 

Our results thus completely classify what can be 
achieved with an n-qubit Hamiltonian and local unitary 
operations. They demonstrate that there are essentially 
only two types of Hamiltonians up to local unitary op- 
erations: those whose Pauli operator expansion contains 
only odd parity terms, and those with at least one even 
term. 

An important caveat to our results concerns efficiency. 
When we state that a set of interactions is universal on a 
set of qubits, we mean simply that those interactions can 
be used to simulate any other interaction, without any 
claim as to whether the simulation procedure is efficient, 
or otherwise. We will say such a set of interactions is dy- 
namically universal. By contrast, in the context of quan- 
tum computing, a set of resources is said to be universal 
for quantum computation if it can be used to simulate a 
standard set of universal gates, such as the controlled- 
not and one-qubit unitary gates, with an overhead that 
is at most polynomial in the number of qubits. 

Our results thus concern dynamic universality, and do 
not directly address the question of universal quantum 
computation. However, some general observations may 
be made about the efficiency of our constructions. When 
the number of qubits that are directly coupled by the 
Hamiltonian is bounded above by some constant k, then 
our simulation techniques for gates such as the controlled- 
NOT only incur an overhead polynomial in the total num- 
ber of qubits. Thus, when the coupling size is bounded, 
our results give criteria not only for dynamic universal- 
ity, but also for universal quantum computation. By con- 
trast, when the number of qubits involved in couplings is 
unbounded, our simulation technique is not polynomial, 
and thus our results cannot be applied to deduce any- 
thing about universal quantum computation. Indeed, we 
conjecture that the two concepts of universality do not, 
in general, coincide. 

For the remainder of this paper, when we speak of a 
set of couplings being universal, we mean dynamically 
universal, unless otherwise stated. The only exception is 
in Section IVII which contains the proof that when the 
terms in the Hamiltonian couple only a bounded number 
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of qubits, our techniques can be made efficient. 

What is the significance of our findings? It is tempting 
to conclude that the main significance is for the design 
of quantum computers. However, we do not believe our 
results are especially significant for such questions. Not 
only are our constructions impractical, but it is an em- 
pirical fact that most interactions occurring in nature 
are two-body interactions, and thus are adequately dealt 
with by earlier work. 

We believe our results are interesting for two other, 
less obvious, reasons. The first is the intrinsic inter- 
est in obtaining general insights into quantum dynam- 
ics. The fact that, given local unitary control, there are 
only two different classes of Hamiltonian evolution, seems 
to us a significant insight into the complicated space of 
possible dynamical evolutions. It tells us that dynam- 
ics within each of these classes are fungible physical re- 
sources. In the case of two-party interactions this insight 
has led to the beginnings of a quantitative theory of the 
strength [l4l Il5| of dynamical operations, much as the 
theory of entanglement dilution and concentration [lq 
led to the quantitative theory of entanglement 0, Il8j|. 

An interesting contrast is to the situation with state 
entanglement, where the multipartite structure is com- 
plex and only partially understood. The number of 
classes of states which are inequivalent under local oper- 
ations and classical communication (LOCC) is enormous 
- for four-qubit states, there are already at least 9 in- 
equivalent classes [Hj • Our results thus demonstrate that 
there are no direct analogies between multiparty state 
entanglement and multiparty entangling dynamics — in 
fact, the situation is substantially simpler for dynamics 
than it is for states. 

A second reason for interest is possible indirect appli- 
cations. For example, although four-qubit interactions 
may not occur in nature, it is certainly the case that 
interactions involving two objects with four-dimensional 
state spaces may occur in nature. Such systems can natu- 
rally be mapped onto our problem by considering a single 
four-dimensional system as being, effectively, a system of 
two qubits. Constructions like this may make our results 
of interest, at least in principle, for realistic physical sys- 
tems. 

The paper is structured as follows. Section [H] pro- 
vides background and definitions. Section ITTT1 establishes 
a body of general techniques for simulating one Hamil- 
tonian with another. These techniques are applied in 
Section llVl to prove our first main result, the character- 
ization of when a Hamiltonian is universal, assisted by 
one-qubit unitaries. Section studies the non- universal 
case, proving that any odd entangling Hamiltonian may 
be used to simulate any other odd Hamiltonian. We then 
describe an encoded universality scheme that allows an 
odd entangling n-qubit Hamiltonian to act universally on 
n — 1 qubits, and provide a Lie-algebraic classification of 
this case. In Section fVll we show that our techniques can 
be made efficient under certain conditions of bounded 
coupling size, and finally in Section IVIII we summarize 



our results. 



II. BACKGROUND AND DEFINITIONS 

This section contains definitions and background ma- 
terial for the remainder of the paper. We begin with some 
notation and definitions, followed by a discussion of what 
it means to simulate one Hamiltonian with another, and 
finally we review some previous work on the universality 
of two-body Hamiltonians. 

We now introduce some notational conventions. As 
stated in the introduction, an arbitrary Hamiltonian H 
on n qubits can be uniquely written in terms of Pauli 
operators via the Pauli operator expansion 
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H = Y2 -3n a h ® • • • ® °j« i (2) 

jl,—,j„=0 

where the hj 1 j n are real numbers and a%, 02, 173 are the 
Pauli sigma matrices X, Y, Z, respectively, with cto = I 
the identity. Let the index a denote each different combi- 
nation ji . . . j n corresponding to non-zero terms in equa- 
tion (J2J. So for each non-zero term in the Pauli expan- 
sion of H we write H a = hj t j <jj 1 ® ... ® Uj and thus 
" V,//„. 

Definition 1. Let C a denote the set of all Pauli product 
Hamiltonians that couple the same set of qubits as H a = 
hj 1 ....j n <Tj 1 ® ... ® crj n , that is C a = {a^ <£> ■■• <£> Cfc„ \h = 
iff ji = 0}. We call each set C a a coupling set. 

For example if H a = X ® / ®Y (g) Y , then C a is the set 
of all products of Paulis acting non-trivially on the same 
qubits, or {X®I®X®X, X®I®X®Y, Z (g) I <g> Z <g> Z} . 

Definition 2. The set S a is the set of qubits coupled by 
H a , or, equivalently, all elements of C a . 

In the example above S a = {1,3,4}. Note that differ- 
ent H a 's can give rise to the same C a and S a . We will 
use the notation \S a \ to denote the number of qubits in 
the set S a . 

Now, what does it mean to simulate one Hamiltonian 
with another? If we can approximately induce evolution 
according to a Hamiltonian H on an arbitrary state \tp) 
for an arbitrary time t without actually using H , then we 
can simulate H , provided the approximation can be made 
arbitrarily good. This concept of simulation is motivated 
by quantum computation which uses a universal set of 
gates to simulate arbitrary unitary evolutions on a set of 
qubits. 

Our approach to the question of whether a Hamilto- 
nian is universal is to exhaustively build up the repertoire 
of different evolutions simulatable with the Hamiltonian, 
in such a way that it becomes clear whether or not the 
repertoire is a universal set of gates. 

A first observation is that, given a term H a = 
hj 1 ...j n o~j 1 ® ■ ■ ■ <£> (Tj n we can simulate xH a for any real 
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non-zero x by adjusting the amount of time we evolve 
according to H a . Therefore, we can always simulate 
H a = ±(Tj 1 ® • • • ® cr Jre with the sign given by the sign of 
hj 1 ...j n . If it is negative, then we can make it positive by 
conjugating H a by a one-qubit unitary U which anticom- 
mutes with it. Thus, we can always obtain aj t ® ■ ■ - ®<Jj n ■ 
For the remainder of the paper, we will always assume 
that terms like H a have this form. 

A simple, but important, observation is that, given the 
ability to evolve according to some Hamiltonian J and to 
perform a unitary operation U and its inverse U' , we can 
evolve according to 

Ue- Ut tf = e~ iUJuh . (3) 

That is, we can simulate evolution according to the 
Hamiltonian J' = UJU'. 

This result can be used to show, for example, that 
given a Pauli product Hamiltonian H a , we can simulate 
any other coupling in C a , simply by performing local 
changes of basis on each of the qubits to interchange the 
role of the x, y and z axes. This is done by conjugating 
by one of the following three rotations: 

e^ x , e^ Y , e^ z . (4) 

Now, suppose we can evolve according to two Hamil- 
tonians Ji and J 2 . Then for small times A the following 
identity holds approximately: 

e -iJiA e — U 2 A ^ e -i(J 1 + J 2 )A /g\ 

That is, we can simulate evolution according to the 
Hamiltonian J\ + J 2 . Equation JSJ) is important because 
it tells us that if we are able to evolve a system according 
to a set of different Hamiltonians, then it is possible to 
simulate arbitrary linear combinations of elements of the 
set. We will treat this identity as though it is exact for 
the remainder of the paper. This is justified for small A. 
(See for an analysis of the errors induced by this ap- 
proximation, and the overhead required to reduce them). 
The above identities are used extensively in this paper, 
as they were in 0, S H S S S S H S 113 

It is also possible to simulate the commutator of two 
Hamiltonians |2Cll2lj. since 

g-tJiAgtJaAgUiAg-iJaA ^ e ~i(i[Ji ,M) A 2 /g\ 

for small A. The error in this approximation is of order 
A 3 , and can be made insignificant by choosing A suffi- 
ciently small. This completes the basic set of tools that 
we use to build up our repertoire of simulatable Hamil- 
tonians. The reason for this is that given a set of Hamil- 
tonians L = { Ji, J^} the set of all simulatable Hamil- 
tonians is given by the Lie algebra generated by the set 
L |22| which can in turn be generated with linear combi- 
nations and i times commutators of elements from L. 

We can thus refine the central question of this paper 
to be: how does the structure of the given Hamiltonian 
H determine the Lie algebra that can be generated by H 
and arbitrary one-qubit evolutions? This is the question 
that we address in the remainder of this paper. 



III. METHODS 

In the previous section we introduced some notation 
and basic tools. This section is concerned with build- 
ing up more sophisticated simulation methods for the 
proofs of our main results, in later sections. In partic- 
ular, there are two interesting simulation ideas — term 
isolation and commutator restriction — that we will ex- 
amine in separate subsections. These ideas may be more 
fully described as follows: 

A. Term isolation: Given that we can simulate H 
and perform arbitrary local unitaries, we can sim- 
ulate an arbitrary term H a in the expansion of H. 
Recall that we write H = J2 a ^ a w here each H a 
is a Pauli product Hamiltonian. 

B. Commutator restrictions: We examine the re- 
strictions placed on the simulation of commutators 
of coupling terms in H . 



A. Term isolation 

In this section, we show that given H we can simu- 
late each term H a in H using one-qubit unitaries and 
the composition identities given in equations @ and JSJ. 
Thus the capacity to simulate H is equivalent up to one- 
qubit unitaries to being able to simulate each H a . 

For simplicity in the proof we now note that we can use 
one-qubit unitaries to simulate the Hamiltonian — 
V\ ® • • • ® V n HVi ® • • • ® V% where the one-qubit unitaries 
V\ , . . . , V n are chosen so that the term H a has all X and 
Y operators in its expansion taken to Z and all Z and / 
operators left alone. Thus in the Hamiltonian every 
<jj 7^ I in H a is now Z . Let us denote this term in 
by H az . This term, equivalent to H a , is the term that we 
wish to isolate. ^From now on we use the convention that 
the superscript on a simulated Hamiltonian indicates a 
step in the algorithm, so would be the Hamiltonian 
simulated after the j th step. We will also use subscripts 
on non-trivial one-qubit operators to indicate which qubit 
they are acting on. For example, X acting on the third 
qubit is written X$ = IiE)I®X®I®---. 

We now use to simulate = Z-lH^Z-l+HW. 
Noting that ZZZ = Z, ZXZ = -X and ZYZ = 
—Y, we see that the term H az is replaced in by 
2H az , but any term in which has X or Y act- 

ing on the first qubit is cancelled out. We then simu- 
late # (3) = Z 2 H^Z 2 + # (2) and so on until we obtain 
H (n+i) = z n H^Z n + ffW, In this Hamiltonian, H az 
has been replaced by 2 n H az and it consists only of terms 
containing Z or /. 

Now we wish to remove all interactions that act on 
qubits outside of S a . We do this by conjugating with X 
on each of these qubits. If H az has an / acting on the q th 
qubit, then we simulate #( n + 2 ) = X q H^ n +^ X q +H {n+1 \ 
This takes 2 n H az to 2 n+1 H az and cancels any terms that 
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have a Z acting on the q th qubit. Let k = \S a \. Since 
there are n — k qubits upon which H az doesn't act, if we 
perform this style of simulation n — k times on different 
qubits, then we will have removed all interactions acting 
on qubits outside of S a . So the simulated Hamiltonian 
jj(2n-k+i) nag j n t erac tions that only act on subsets of 
S a , and H az has become 2 2n ~ k H az . 

At this point, we wish to eliminate all remaining terms 
except H az . Denote one of these undesirable terms 
by Hp z . We know that Hp z couples a set of qubits 
5^3 that is strictly contained in S a . Thus, there must 
be some qubit q which is in S a and not in Sp. Not- 
ing that XpXq commutes with Z p Z q but anticommutes 
with Zpl q for p ^ q, we see that conjugating by X p X q 
leaves H az invariant but takes Hp z to —Hp z if we choose 
p to be in Sp. Thus, if we simulate i/( 2 ™- fc + 2 ) = 
X p X q H^ n - k+1 1X p X q +H^ n - k+1 ^ we eliminate H 0Z . It- 
erating this procedure for every combination of p and q 
in S a , we are can eliminate every remaining undesirable 
term. There are ( \ ) such possible combinations, so we 
finally obtain 

fl .(2n-fc+l+( \ )) = 2 2n-fc+( \ ) R ^ ^ 

Summarizing, we can isolate any term in H, and thus 
can simulate all elements of C a for every H a appearing 
in H, as well as all linear combinations of the elements 
of C a . 

We note in passing a group-theoretic interpreta- 
tion j2^, 13 °f the term isolation procedure described 
above. A particular term H a that we wish to isolate 
forms, along with the identity, an order two group Q 
which is a subgroup of the full Pauli group V on our 
system. In particular H a and I are a representation of 
this subgroup Q . Denote the commutant of Q in V (the 
set of all group elements in V which commute with Q) as 
Q 1 . If we map the elements of Q 1 to Pauli operators, then 
we have a faithful representation of Q' . Let us denote the 
elements of this representation by D(g'), where g' € Q. 
From Schur's lemma it then follows that by averaging 
over all elements in this representation of Q' we obtain 
only elements in the representation of Q given by H a and 
I: 

T^TT £ D(g')HD(g'y = al + bH a (8) 

for some constants a and b. Thus we can isolate a term 
by performing the appropriate group average over the 
commutant subgroup. In fact, since Q is abelian, we can 
average over elements in Q' which are not in Q. The 
term isolation procedure we described above is a concrete 
realization of this group average. 

B. Commutator restrictions 

We now turn to the simulation of Hamiltonians using 
commutators, focusing on the possible forms of Hamil- 



tonians simulated in this manner. The restrictions we 
obtain will be vital for the results of the next section. 

Consider two different elements of the coupling set C a , 
H a and H a i. It is straightforward to verify from the 
commutation relations for the Pauli matrices that the 
commutator [H ai H a i] is non-zero if and only if there is 
an odd number of locations in S a where the couplings H a 
and H a i differ. For example, consider the commutator 
of H a = X ® X (g) X and H a , = Y <g> X ®Y, both of 
which couple the same qubits: 

i[H a ,H a >] =i(XY <g> I <g> XY -YX <g> I ®YX) 

=i(XY ®I®XY-XY®I® XY) = 0. 

(9) 

We see that in this case an even number — two — of 
the qubits are acted on by different Paulis, and so the 
commutator is zero. However if H a i = Y ® X ® X, 
the commutator is non-zero, indeed it is i[H a , H a i] = 
— 2Z ® I ® I. We see from this example that given an 
initial coupling, we can generate terms which couple a 
different set of qubits. 

Suppose now that C' a and Cp are two different cou- 
pling sets. If they act on nonoverlapping sets of qubits, 
then any commutator between an element of C a and an 
element of Cp will always be zero. The commutators be- 
tween two particular Hamiltonians H a and Hp depend 
only on their actions on qubits in the intersection of the 
sets S a and Sp. Recall that the commutator of two terms 
is nonzero if there is an odd number of pairs that dis- 
agree. Thus H a and Hp must differ on an odd number 
of qubits from the intersection of S a and Sp in order to 
have a nonzero commutator. To summarize, a commu- 
tator of Hamiltonians from C a and Cp can generate a 
Hamiltonian coupling a set of qubits S-y precisely when 
S-y is in S a U Sp, S 7 contains an odd number of qubits 
from S a ("I Sp, and S-y contains all the qubits from both 
S a and Sp which are not in S a ft Sp. 

This becomes clearer with an example. Suppose H a = 
X X ® X (g) X and Hp = Z®X®X®I. We find 
that the commutator is i[H a , Hp] = —2Y 
So we see that we have simulated a two-qubit entangling 
Hamiltonian, with a four- and a three-qubit coupling. 

Combining all of our results about simulation so far, 
we see that given H — ^ Q H a and arbitrary one-qubit 
unitaries, we can isolate any term H a . This can then 
be used to simulate any Hamiltonian coupling the same 
set of qubits. All of these terms can be combined in ar- 
bitrary linear combinations. Finally, we can use pairs of 
couplings C a and Cp to simulate a coupling on a different 
set of qubits 5 7 when the conditions stated above hold. 
The coupling C 7 can be added to our repertoire of sim- 
ulatable couplings, and can be used in turn to generate 
new couplings. How many different couplings are there? 
Since we assume that there is a finite number of qubits n 
there are no more than 2™ different coupling sets, so the 
process of generating new couplings and adding them to 
the repertoire must terminate. 
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Once all of the simulatable couplings have been enu- 
merated, the complete set of simulatable Hamiltonians 
consists of those whose Pauli operator expansion con- 
tains only terms that belong to one of the simulatable 
couplings. Of course, this procedure of exhaustive enu- 
meration for determining which Hamiltonians are simu- 
latable is not especially efficient or insightful. In the next 
section we provide a surprisingly simple procedure that 
enables us to determine when a Hamiltonian is universal. 



IV. STRUCTURE OF SIMULATED 
COUPLINGS AND DYNAMIC UNIVERSALITY 

In the previous section we identified various simula- 
tion methods for Hamiltonians acting on qubits. In this 
section we use these methods to classify which many- 
qubit Hamiltonians are universal given local unitary op- 
erations, assuming throughout that the Hamiltonians un- 
der consideration are entangling. In order to find this 
classification we use properties of the Pauli operator ex- 
pansion of a qubit Hamiltonian as given by equation J2J . 
In particular, we will see that the parities of the cou- 
plings in this expansion determine whether or not H is 
universal. 

We know from subsection llll Al that given H — J2 a H a 
and local unitaries we can simulate any particular cou- 
pling term H a . If we had as our base set of operations 
each H a in the expansion of H, and local unitaries, we 
could simulate H. Thus, having H and local unitaries 
is equivalent to having {H a } and local unitaries. In this 
section we focus on the simulating capacity of particular 
coupling terms H a as this will be sufficient for analyzing 
the universality of H . 

In section [HI we defined the coupling set C a as be- 
ing the complete set of Pauli product Hamiltonians that 
couple the set of qubits S a . We also noted that a single 
element of C' a and arbitrary one-qubit unitary control 
generates all elements of C a . In this section we will often 
use the coupling set C a and any single element of that set 
interchangeably. For convenience, we write LU to rep- 
resent all local unitaries, that is, products of one-qubit 
unitaries. We also define the parity of a Pauli-product 
Hamiltonian H a to be odd if it acts non-trivially on an 
odd number of qubits, or, equivalently, if S a contains an 
odd number qubits. Otherwise, we say that H a has even 
parity. We will see that dynamic universality of H is 
completely determined by the parities of the terms in H . 

For convenience, we restate a previous result that we 
will use frequently 

Theorem 1 (Bipartite Hamiltonian theorem [l], 0, 

0, 0, Ht ) • Suppose H , acting on n qubits, has only one- 
and two- qubit terms in its Pauli-product expansion, and 
that H is entangling, i.e., all n qubits are connected, pos- 
sibly indirectly, by the terms in H . Then H , together with 
local unitaries, is universal for quantum computation on 
n qubits. 



We now prove a series of lemmas leading to our first 
new theorem. 

Lemma 1. // Hamiltonians H a and Hp act on sets of 
qubits S a and Sp such that Sp C S a and \S a \ = \Sp\ + 1 
then the set {H a , Hp, LU} is universal on S a . 

Before giving the proof, consider the example H a = 
X ® X X (g) X and Hp = Y X X ® I. Their 
commutator is 

i[H a ,Hp]=-2Z®I®I®X (10) 

which means that we can simulate arbitrary couplings 
between qubits 1 and 4. On the other hand, the commu- 
tator of X® Y®X®I € Cp with H a generates a coupling 
between qubits 2 and 4. Similarly, we can couple qubits 
3 and 4. Using Theorem \I\we see that {H a , Hp, LU} is 
universal on the set of qubits S a . 

Proof: Assume without loss of generality that the 
qubits are numbered so that the first n — 1 of them are 
in Sp. Then X® n G C a and X® n ~ 2 ®Y ® I £ Cp, 
where Y acts on the (n — l) th qubit. We can simulate 
the commutator of these Hamiltonians 

i[X® n , X® n ~ 2 <g> Y <g> I] = -21® n - 2 ®Z®X. (11) 

and thus we can couple the (n — l) th and nth qubits. We 
can perform similar simulations where Y acts on each 
qubit in the range 1 to n — 1. This generates two-qubit 
couplings connecting all of S a , and Theorem ^ implies 
that {H a , Hp, LU} is universal on S a . □ 

Lemma 2. If a Hamiltonian H a has even parity then 
{H a ,LU} is universal on S a . 

Proof: This result follows almost immediately from 
Lemma ^ Let n = \S a \. Notationally, it will be conve- 
nient to omit qubits outside the set S a in the following; in 
all cases there is an implied identity action on the omit- 
ted qubits. Since X® n G C a and y®"" 1 ® X & C a , we 
can simulate the commutator 

i[X® n ,Y® n - 1 ®X] = 2i{iZ)® n ~ 1 ®I (12) 

as n — 1 is an odd number. Since 2i(iZ)® n ~ 1 ® / acts 
on n — 1 qubits, Lemma 1 allows us to conclude that 
{H a , LU} is universal on S a . □ 

Corollary 1. // H a and Hp both have even parity and 
S a nSp 7^ 0, then {H a , Hp, LU} is universal on S a USp. 

Proof: Lemma |21 tells us that H a is universal on S a 
and Hp is universal on Sp. Since they have a non-empty 
intersection, they can simulate a set of two-qubit Hamil- 
tonians connecting every qubit in S a U Sp which implies 
that {H a , Hp,LU} is universal on 5 a U 5,3. □ 

Lemma 3. If H a has odd parity, n, Hp has even parity, 
and Sp C S a , then {H a , Hp, LU} is universal on S a . 
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Consider the following example. Let H a — X® 5 and 
Hp = J® 3 The commutator of these Hamilto- 

nians is a four-qubit coupling, and therefore, by Lemma 
^ they are universal on S a : 

i[H a ,Hp] = -2X®X®X®Z®I (13) 

Let's now generalize this example to prove Lemma [3] 

Proof: Using H a and LU we can simulate the Hamil- 
tonian X® n , and by Lemma^Hp and LU can be used to 
simulate the Hamiltonian I®" -2 ® Y ® (assuming we 
number the qubits so that the n th and (n — l) th qubits 
are in S a ). Then we can simulate the Hamiltonian 

i[X® n , I® n ~ 2 ®Y®X] = -2X® n ~ 2 ®Z®I, (14) 

which acts on n — 1 qubits in S a . Therefore, Lemma Q 
implies that {H a , Hp, LU} is universal on 5„. □ 

Corollary 2. If H a has odd parity and we have a uni- 
versal set of gates Up acting on a set of qubits Sp such 
that Sp C S a and \Sp\ > 1, then the set {H a , Up, LU} is 
universal on S a . 

A simple example of such a set of universal gates on 
three qubits is {X m , CNOT <g> I, LU}. 

Proof: If we have a universal set of gates Up on Sp, 
then it is possible to simulate a unitary operator equiva- 
lent to a Hamiltonian evolution by H 7 acting on an even 
number of qubits \Sj\ with 1 < \Sj\ < \Sp\. Then by 
Lemma |3 the corollary is true. □ 

Lemma 4. If H a has even parity, n, Hp has odd parity, 
m, and S a C\Sp ^ 0, then {H a , Hp,LU} is universal on 
S a U Sp. 

Proof: First, consider the case where \S a H Sp\ = 1. 
Label the qubits so that S a contains the first n qubits 
from the left and Sp contains the first m qubits from the 
right. Thus x®"®/®" 1 - 1 e C a and j®"-i<g,y®™ g Cp. 
We can simulate the commutator 

Now, this commutator acts on |S' Q US , / 3| = m+n—1 qubits 
which is an even number, and therefore by Lemma|5]this 
Hamiltonian and LU are universal on these qubits. 

The case \S a fl Sp\ > 1 is even simpler to prove. Since 
n is even, H a is universal on S a by Lemma|21 Now, when 
\S a H Sp\ > 1 we have a universal set of gates acting on 
S a C\Sp, so Corollary|21proves that we have a universal set 
on Sp , and therefore we have a universal set on S a U Sp. 
□ 

Lemma 5. If we have a universal set of gates U a acting 
on S a such that \S a \ > 2, and Hp acting on Sp such that 
S a (~]Sp ^ 0, then {U a , Hp, LU} is universal on S a USp. 

Proof: The proof follows simply by considering the 
four possible parity combinations for \S a \ and 

• Case: \S a \ even, \Sp\ odd. The result follows from 
Lemma 0] 



• Case: \S a \ even, \Sp\ even. The result follows from 
Corollary ^ 

• Case: \S a \ odd, \Sp\ even. The result follows from 
Lemma 0J with the roles of a and (3 interchanged. 

• Case: \S a \ odd, \Sp\ odd. Pick an even parity 
subset of S a which has a non-trivial overlap 
with Sp. Since we have universality on S a we must 
also have universality on S ev . Lemma 0] therefore 
implies universality on S ov U Sp. Universality on 
S a U Sp now follows from Theorem ^ □ 

Lemma |S] gives us a composition rule for determining 
whether or not a Hamiltonian is universal on a given set 
of qubits. It tells us that if we have a set qubits cou- 
pled by even-parity coupling terms (and thus a universal 
set of gates on the set), then any other qubits that they 
are connected to, even indirectly, will also have a univer- 
sal set defined on them. We will see from the following 
lemma that this is the only way a universal set of gates 
can be derived. 

Lemma 6. If H a and Hp both have odd parity, then 
their commutator H 1 = i[H a ,Hp] is either or it has 
odd parity. 

Proof: If S a n Sp = then i[H a ,Hp] = 0, and the 
lemma is trivially true. In the case where S a DSp ^ let 
u = \S a nSp\. Now, Hj ^ only when |S , 7 n(S , ct nS , /3 )| = d 
is odd. By definition we find that \Sj\ = \S a \ + \Sp\ — 
2u + d. Since \S a \ + \Sp\ and 2u are even, but d is odd, 
|5 7 | must be odd for any non-zero H~ f . □ 

The conclusions of these lemmas can be succinctly ex- 
pressed in the following theorem. 

Theorem 2. An entangling Hamiltonian acting on n 
qubits, together with local unitary operations, is dynami- 
cally universal if and only if the Pauli operator expansion 
of the Hamiltonian contains a term with an even number 
of entries. 

Proof: From Lemma |S] we know that an entangling 
Hamiltonian whose Pauli operator expansion contains 
only odd-parity terms, together with local unitaries, 
can only generate other odd-parity Hamiltonians, which 
shows that it is not universal. Conversely, from Lem- 
mas |2 and El we see immediately that an entangling 
Hamiltonian with at least one even-parity term is uni- 
versal. □ 

We now know that in order for a general n-qubit 
Hamiltonian to be universal when aided by one-qubit uni- 
taries, it must fulfill two conditions. The first condition 
is that the Hamiltonian must be entangling, in the sense 
explained in the introduction, i.e., all n qubits must be 
connected, either directly or indirectly, by coupling terms 
in the Hamiltonian. The second condition is that at least 
one of the coupling terms must have even parity. Thus, 
our result reduces the problem of determining when a 
Hamiltonian is universal to that of counting the parity of 
terms in its Pauli operator expansion. 
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Let's look at a couple of examples. It was suggested 
in the preprint version of pj (but not the published 
version) that it might be possible to construct many- 
body Hamiltonians which are not universal by using re- 
sults from the theory of entanglement. For example, 
the GHZ state \GHZ) = (jOOO) + |lll))/\/2 and VP- 
state (|001) + |010) + |100))/V3 are distinct types of 
entanglement which cannot be interconverted pQj. even 
stochastically, by local operations and classical com- 
munication. This led Q] to conjecture that Hamilto- 
nians such as Hghz = \GHZ)(GHZ\ and Hqhz' = 
\GHZ) (000| + 1 000) (GHZ\ are not universal, when as- 
sisted by local unitaries. This conjecture turns out to be 
incorrect. Expanding in the Pauli basis and omitting (g) 
for brevity, we obtain: 

H GH z oc III + ZZI + ZIZ + IZZ — XYY 

-YXY - YYX (16) 
H GH z' oc III + ZII + IZI + II Z + ZZI 

+ZIZ + IZZ + ZZZ - XXX. (17) 

In both cases, we simply check that each qubit is coupled 
by one of the two- or three-qubit terms, so the Hamilto- 
nian is entangling, and note that there are terms with 
even parity, so by Theorem [5] both these Hamiltonians 
are universal when assisted by local unitaries. 



V. HAMILTONIANS WITH ALL-ODD PARITY 

In the previous section we showed that the only non- 
universal entangling Hamiltonians are the odd Hamilto- 
nians, i.e., those whose couplings all act on an odd num- 
ber of qubits. In this section we study what dynamical 
operations can be achieved using such Hamiltonians, to- 
gether with local unitary operations. We prove two main 
results. 

The first result is that an odd entangling Hamiltonian 
can be used to simulate any other odd Hamiltonian on 
the system of n qubits. This result, in combination with 
the results of the previous section, shows that there are 
essentially only two types of entangling Hamiltonian on 
n qubits: Hamiltonians that are odd, and those that are 
not. Within these two classes all the Hamiltonians are 
essentially inter-convertible, in the sense that any one 
can be used to simulate the other. Furthermore, the en- 
tangling Hamiltonians that are not odd are intrinsically 
more powerful than the odd Hamiltonians, since they can 
be used to simulate any odd Hamiltonian, but not vice 
versa. 

The second result is to show that odd Hamiltonians 
can be made universal, by using an appropriate logical 
basis of qubits for our system, similar to the ideas used 
in quantum error-correction. In particular, we show that 
such an interaction on n qubits is universal on a set of 
ra—1 logical qubits. In fact, we will see that the encoding 
is as simple as it could be: each of the n — 1 logical qubits 



corresponds directly to one of the original n qubits, while 
the single qubit left over is not used. 

To obtain our results we first need a simple lemma 
allowing us to use an odd Hamiltonian coupling a set of 
qubits to generate odd Hamiltonians acting on a subset 
of those qubits. 

Lemma 7. If H a has odd parity then we can simulate 
any other odd-parity Hamiltonian Hp provided Sp C S a . 

Proof: Let \S a \ = 2m +1. We prove this lemma using 
induction on m. The lemma is trivially true for the case 
m = 1 (that is, \S a \ — 3). Now, we make the inductive 
assumption that the lemma holds for the m th case and 
prove that it holds for the (m -I- l) th case. Given H a i 
acting on a set of \S a '\ = 2m + 3 qubits, we need only 
show that we can simulate a Hamiltonian acting on any 
subset of S a > containing 2m + 1 qubits. By assumption, 
we can simulate the (2m+3)-qubit Hamiltonians X® 2m+3 
and X® 2 ® Y 2m+1 and thus their commutator 

z[l® 2m+3 ,l 82 ®y 2m+1 ] = -2iI®I®{iZ)® 2m+1 . (18) 

Therefore we can simulate a Hamiltonian acting on the 
final 2m + 1 qubits of S a . Similarly, we can simulate a 
Hamiltonian acting on any subset of S a containing 2m+l 
qubits, which proves the inductive hypothesis and thus 
the lemma. □ 

It turns out that both the main results of this sec- 
tion are corollaries of this lemma and a second result 
that is motivated by the following example. Suppose 
we have the ability to simulate the odd Hamiltonian 
H = ZZZII+IIZZZ (where we have omitted both sub- 
scripts and tensor products). We already know that we 
can therefore simulate ZZZII and IIZZZ separately. 
Suppose further that we wish to simulate another odd 
Hamiltonian, ZIIZZ. We can see that this is possible, 
using our intuition about non-odd Hamiltonians, by "iso- 
lating" one of the qubits, say the fifth qubit, and consid- 
ering the Hamiltonians that we can simulate on the first 
four qubits. So, let's alter our example and give ourselves 
the ability to to perform ZZZI and IIZZ on the first 
four qubits, and we will attempt to simulate ZIIZ. No- 
tice that now we have an odd and an even term that con- 
nect the first four qubits, so by Theorem |2 there must 
be a sequence of commutators and linear combinations 
that allow us to simulate ZIIZ. Here is such a sequence, 
where in each step we generate a new coupling to add to 
our set of allowed couplings: 

i[ZZYI,IIXZ] = 2ZZZZ, 
i[XZXX,YZYY] = 2ZIZZ, (19) 
i[ZIZY,IIZX] = 2ZIIZ. 

Now, if we consider the original problem on five qubits, 
we see that the odd parity restriction shows that this 
procedure generates the desired coupling ZIIZZ: 

i[ZZYII,IIXZZ] = 2ZZZZZ, 
i{XZXXZ, YZYYZ] = 2ZIZZI, (20) 
i[ZIZYI,IIZXZ] = 2ZIIZZ. 
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In a similar vein, suppose that we wish to simulate 
an even-parity coupling ZIIZI. We know that this is 
not possible directly (otherwise H would be dynamically 
universal), but it is possible using a very simple encod- 
ing. If we place the fifth qubit in the Z-eigenstate |0), 
then our procedure above for simulating ZIIZZ allows 
us to simulate ZIIZ on the first four qubits. Thus, H is 
dynamically universal on the first four qubits. 

Our main results of this section generalize these two 
examples. Both examples rely crucially on the fact that 
there was a qubit that could be isolated, that is, acted on 
by only a single member of the set of couplings that we 
used to do our simulation. For example, suppose we tried 
to use the same approach to show that H can simulate 
ZIIZ on four of the five qubits by placing the third qubit 
in the state |0). Then the two couplings that we have at 
our disposal on the remaining four qubits are ZZII and 
IIZZ, which do not connect the four qubits, and thus 
cannot be universal on them. We formalize this intuition 
in the following lemma. 

Lemma 8. Let H = Y] H a be an odd entangling n-qubit 
Hamiltonian, and TL = {H a } be the set of all terms in H . 
Then there exists a set Ai C H that (a) connects all n 
qubits in such a way that (b) at least one of the qubits is 
only acted on by a single element of Ai . We call such a 
qubit an isolated qubit, and Ai an isolating set for that 
qubit. 

Proof: We prove this lemma by giving a constructive 
procedure to generate an appropriate set Ai. For conve- 
nience, define n a = \S a \ for all a. 

1. Choose a term Hi from the set TL. Without loss 
of generality, we may number the qubits so that it 
acts on the first n\ qubits. Add H\ to the set Ai. 

2. Search for a second term in the set TL that overlaps 
with Hi and also acts on at least one qubit outside 
of Si. If there is no such term, then Hi must couple 
all of the qubits, in which case Ai — {Hi} satisfies 
the conditions above and we are done. 

3. Otherwise, choose such a term and call it Hy . De- 
fine ns 1 ns 2 , to be the number of qubits in Si flSV- 
Without loss of generality, we may assume that 
these qubits are strung out in a line, with the ni 
left-most qubits in Si and the right-most qubits 
in 52', and the ns in s 2 , overlapping qubits in the 
middle. 

4. If ns 1 ns 2 i is odd, then use Hy to simulate a Hamil- 
tonian Hi that acts on qubits rc.i,...,ni + — 
n SinS 2 , (where we number from the left, starting at 
1). i?2 acts on the right-most ni = — ns 1 ns 2 , + 1 
qubits, overlapping with Hi on just a single qubit 
(the nith qubit). This is possible since is odd 
and S2 C S2 1 , satisfying the conditions of Lemma0 
Add #2 to M 



5. On the other hand, if ns in s 2 , is even, then we use 
i?2' to simulate a Hamiltonian H2 that acts on 
qubits ni — 1, ...,m + — n-s^S,,- This H2 acts 
on the right-most n% = ri2> — n>Sir\S 3 i + 2 qubits, 
overlapping with Hi on just two qubits, in posi- 
tions rii — 1 and m. Again, 712 is odd and S2 C S2 1 , 
satisfying the conditions of Lemma [3 Add H2 to 
M. 

6. Now, if there are no other terms that overlap H2, 
then the right-most qubit in S2 must be isolated 
since H2 acts on at least three qubits, and overlaps 
with Hi on at most two. If we then add the remain- 
ing Hamiltonians from TL (i.e., all except for Hi and 
H2') to Ai, then Ai must couple all n qubits and 
contain an isolated qubit, satisfying the conditions 
above, and so we are done. 

7. Otherwise, repeat steps 2 to 5 to generate a Hamil- 
tonian H3 that overlaps only with H2 on one or 
two qubits, and add it to Ai. Repeat this process, 
adding a Hamiltonian to Ai each time, until it be- 
comes impossible to find a term that both overlaps 
with the previous term and acts on at least one 
more qubit than it. When the process terminates 
(which must happen eventually since there is only 
a finite number of qubits) , the last term that was 
added must contain an isolated qubit. If Ai con- 
nects all n qubits then we are done, otherwise add 
the remaining Hamiltonians from TL to Ai to com- 
plete the construction. 

□ 

Using this lemma, we can prove our two main results. 

Theorem 3. Let H be an odd n-qubit entangling Hamil- 
tonian. Then H and LU can simulate any odd Hamilto- 
nian on the n qubits. 

Proof: Let Ai be an isolating set for a qubit, which, 
without loss of generality, we may choose to be the n th 
qubit. Suppose we consider the set of couplings Ai' on 
the first n — 1 qubits that arises by simply taking the 
couplings in Ai, and omitting the Pauli acting on the 
final qubit. Then by construction of Ai we see that this 
set (a) connects the first n — 1 qubits, and (b) contains 
an element that acts on an even number of these qubits, 
corresponding to the element of Ai that couples to the 
isolated qubit. By Theorem[21it follows that Ai', together 
with local unitaries, is universal on the first n — 1 qubits. 

Lifting back up to the full set of n qubits, we see that 
Ai must generate the set of all odd couplings on the n 
qubits. To see this a little more explicitly, suppose that 
we wish to generate an odd coupling er. Let a' be the 
corresponding coupling on the first n — 1 qubits. By an 
appropriate sequence of commutators of elements of Ai' 
we can generate a . The corresponding sequence of com- 
mutators in Ai must generate a, up to possible relabeling 
on the final qubit, which can be accomplished via appro- 
priate local unitaries. □ 
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Theorem 4. Let H be an odd n-qubit entangling Hamil- 
tonian. Then H and LU are universal on a set of n — 1 
logical qubits. 

Before proving this theorem, let's consider another ex- 
ample. Consider the Hamiltonian H a = Z ® Z ® Z . We 
know from the previous section that this Hamiltonian 
and local unitaries do not form a universal set of opera- 
tions on three qubits. However, if we prepare the third 
qubit in the |0) eigenstate of Z at the beginning of the 
procedure, then any succession of evolutions of H a and 
local unitaries acting only on the first two qubits will 
leave the third qubit invariant throughout the evolution. 
So if, for instance, we trace over the third qubit we find 
that the reduced system evolves according to Z ® Z. We 
know that this Hamiltonian is universal on the reduced 
system, so in effect we have a universal set of gates on the 
first two qubits. Essential to this example is our ability to 
identify a qubit that can be prepared in a local eigenstate 
of the Hamiltonian that entangles the other qubits. 

Proof: The proof is simply to number the qubits 
1, . . . ,n, and to prepare the n th qubit in a fixed state 
|0). Suppose now that we wish to simulate an arbitrary 
Pauli a acting on the first n — 1 qubits. If a has odd par- 
ity, then we can use the results of Theorem|3|to simulate 
a directly. If a has even parity then we use the results 
of Theorem |3 to simulate er ® Z, which leaves the state 
of the final qubit unchanged, and evolves the first n — 1 
qubits according to the Hamiltonian a. Thus, we can use 
this construction to simulate any interaction on the first 
n — 1 qubits. □ 

We now understand that odd n-qubit Hamiltonians, in 
contrast to the non-odd Hamiltonians, do not generate 
the algebra su(2"). What algebra do they generate? 1 It 
turns out that the answer depends on whether n is odd 
or even. 

To state and prove our results, it is helpful to be a little 
more precise about the various Lie algebras we are con- 
sidering. We define g to be the real Lie algebra generated 
by odd parity Paulis acting on n qubits. More precisely, 
2 is a real vector space whose basis elements are of the 
form ia, where a is an odd parity Pauli. We have shown 
that this is the relevant Lie algebra associated with an 
entangling Hamiltonian on n qubits, plus one-qubit uni- 
taries. The following theorem relates g to the standard 
classification of Lie algebras: 

Theorem 5. Ifn is even, then q is isomorphic toso(2 n ). 
If n is odd, then g is isomorphic to sp(2 n ). Furthermore 
the representation of 2 provided by tensor products of odd 
parity Paulis on n qubits is the fundamental (defining) 
representation of these Lie algebras. 

Amusingly, the n = 2 case of this theorem is a well- 
known result from Lie theory, the isomorphism between 



so(4) and su(2)®su(2). This result has received wide use 
in quantum information theory in a different guise — it 
is just the fact that local (special) unitary operations on 
two qubits correspond to real orthogonal transformations 
in the so-called "magic basis" . 

Proof: We consider the n even case first. Let us define 
an operation f(a) = (— l) wt ' CT '<7 on Pauli matrices, where 
wt(cx) is the weight of a. This operation can be extended 
by linearity to all matrices. Observe that: 

f{A) =Y® n A T Y® n . (21) 

Then the Lie algebra q consists of all matrices A such 
that: 

f(A) = -A, and A* = -A. (22) 

The Lie algebra so(2 n ) can be defined similarly. Recall 
that the defining representation of the Lie algebra so (2™) 
consists of 2" x 2" matrices B which satisfy [2^, [2(J 

B T = —B, and B f = —B, (23) 

where T denotes the transpose operation. 

We aim to find a unitary U such that B satisfies 
Eq. (ED if and only if A = UBU f satisfies Eq. (j22J. 
It is easy to see that for any U, A^ — — A if and only if 
B^ = —B, so we need only find a U such that f{A) = —A 
if and only if B T — —B. Straightforward algebraic ma- 
nipulation shows that U = (I — iY® n )/^/2 satisfies this 
requirement. 

The n odd case is very similar. The defining repre- 
sentation of the Lie algebra sp (2™) consists of 2™ x 2" 
matrices which satisfy |25l [2(j 

J^B T J = —B, and £ f = —B, (24) 

where J = Y\ ® I is the Pauli Y acting on the first qubit 
alone. Setting U = I (I - iY®™- 1 )/^, we see that 
Eq. is equivalent to Eq. (J23) if we set A = UBW. 
□ 

An interesting consequence of Theorem [5] occurs for 
odd parity Paulis acting on an odd number of qubits, 
n. Suppose we are given some known pure state, 
We can then ask the question of whether we can trans- 
form this state into any other state, \cf>). Clearly if we 
have control over su(2") we can perform this task. A 
theorem from the study of the controllability of quantum 
systems [27], shows that this task can be performed 
for all states l^) and \<p) if and only if one has control over 
the Lie algebra su(2") or the Lie algebra sp(2"). Thus, 
while we do not have full unitary control when we have 
odd parity Paulis acting on an odd number of qubits (ex- 
cept when n — 1), we can transform any state into any 
other state using these operations. 

VI. EFFICIENCY 



This question was originally posed to us by Greg Kuperberg. 



We have examined the problem of Hamiltonian simu- 
lation across sets of qubits that are coupled by a fixed 
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natural Hamiltonian. However, our results appear to be 
limited by the fact that our simulation techniques are 
manifestly exponential in the total number of qubits n. 
The problem is that our procedure for isolating a single 
term of an arbitrary Hamiltonian requires on the order 
of 2 n local unitaries to be interspersed in the evolution at 
different times (see equation 100 • This is in sharp con- 
trast to the case of Hamiltonians whose Pauli expansions 
contain only two-qudit couplings, and simulation tech- 
niques are polynomial in the number of qudits. This de- 
fect could be fixed if we were to find efficient techniques 
for term isolation, for all other steps in our procedure 
were efficient, at least in principle. 

It is not surprising that, in this very general situ- 
ation, we have not obtained efficient simulation tech- 
niques. Suppose, for example, that we have a family 
of Hamiltonians that has the property that the member 
of the family that acts on n qubits has a decomposition 
containing a tensor product of Xs acting on every subset 
of those n qubits: 

H{2) = XiX 2 

i?(3) = X1X2 + X1X3 + X2X3 + X1X2X3 

i?(4) = X1X2 + ■ • • + X1X2X3 + ■ • • + X1X2X3X4 

(25) 

where the subscripts in brackets indicate the number if 
qubits acted on by each Hamiltonian. These Hamiltoni- 
ans are not at all natural since every possible coupling 
of qubits is represented, regardless of how "far away" the 
qubits are from one another. With such a large number of 
couplings (2" — n — 1 for the n-qubit member of the fam- 
ily), it is not surprising that we have not found a method 
to efficiently simulate a set of gates that is universal for 
quantum computation, since it is difficult to "turn off" 
enough unwanted interactions, while retaining computa- 
tional universality, even though these Hamiltonians are 
dynamically universal. As mentioned in the introduction, 
we conjecture that a generic entangling Hamiltonians will 
not be universal for quantum computation. 

On the other hand, it would be surprising if the follow- 
ing family of Hamiltonians could not be used simulate a 
universal set of gates: 

H(2) = X1X2 

-0(3) = X1X2 + X2X3 + X1X2X3 

-"(4) = X1X2 + X2X3 + A3A4 + X\ X2X3 + A2A3A4 

H(5) = XiX 2 + • • • + X4X5 + X1X2X3 + • • • + A3A4A5 

(26) 

In these Hamiltonians, each qubit is directly coupled to 
at most four other qubits by terms acting on at most 
three qubits. For example, qubit 3 only ever couples di- 
rectly to qubits 1, 2, 4, and 5. Furthermore, the number 



of coupling terms is linear in the number of qubits - 
the Hamiltonian on n qubits has only 2n — 3 terms in its 
decomposition. This number is sufficiently small, by con- 
trast with the general case, that we might hope that it is 
possible to turn off most (but not all) of these couplings 
in an efficient fashion. We call a Hamiltonian (or, more 
precisely, a family of Hamiltonians) k-local if the maxi- 
mum number of qubits coupled by any term in its decom- 
position is k and if the absolute values of the non-zero 
coefficients multiplying each coupling in H (the hj 1 j n 
in the expansion of H in equation J2J) are bounded be- 
low by a constant. This family of Hamiltonians is thus 
3-local. 

Motivated by these examples, we now describe a ran- 
domized algorithm which shows that a single term in a k- 
local Hamiltonian can be isolated with high probability in 
a number of steps that is exponential in k but polynomial 
in n. Thus, given a family of fc-local Hamiltonians, for 
some fixed k, we have a procedure to efficiently simulate 
a set of gates that is universal for quantum computation. 

More precisely, suppose we wish to isolate a single term 
H a in the expansion of H using local unitaries. We give 
a procedure to use a randomly chosen set of local uni- 
taries {Uj} to isolate H a with failure probability bounded 
above by N/2 m , where N is an upper bound on the num- 
ber of terms in H and m is the number of local unitaries 
in the set {Uj}. 

If H is fc-local, then N is polynomial in n — a simple 
bound on N is n k — and hence the probability of failure 
is polynomial in n and decreases exponentially in the 
number of local unitaries. More precisely, to bound the 
probability of failure to be less than e, it turns out to 
be sufficient to choose m > \og(N/e). The number of 
timesteps required increases by a factor of two for each 
extra unitary, so the number of timesteps is bounded 
above by 2 m = N/e. Provided N is a polynomial function 
of n, as is the case if H is fc-local, the number of timesteps 
is also polynomial in 11. 

Algorithm: To explain the algorithm, we begin by 
explaining how to eliminate a single unwanted term from 
H, without worrying about keeping our desired term. 
(We will see later that a simple modification of this pro- 
cedure eliminates the unwanted term and keeps the de- 
sired term.) Suppose, without loss of generality, that the 
unwanted term has the form Hp = Z\Z 2 ■ ■ ■ ZiIi + \ ■ ■ ■ I n . 
In order to eliminate this term, we choose a set of local 
unitaries {Uj}, each of which is a tensor product of n 
unitaries each chosen independently with equal probabil- 
ity from the set {I,X,Y, Z}. We then do the following 
conjugations: 

H (1) = U x HUi + H, 

= u 2 hWu 2 + hW, 

(27) 

= U 3 H {m -^Uj + ff( m - x ) 
The probability that £f( m ) still contains the unwanted 
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term Hp is the probability that each Uj commutes with 
Hp . The probability that a particular Uj commutes with 
Hp is simply the probability that the total number of X 
and Y terms in Uj that act on the first I qubits is even. 
So, if I is 1, then the probability that Uj commutes with 
Hp = Z\ is just the probability that Uj has a Z or an 
I on the first qubit, which is 1/2. Similarly, if Hp acts 
nontrivially on the first 2 qubits, then the probability 
that Uj commutes with it is the probability that its first 
two terms are I\I 2 , Z X Z 2 , X X X 2 , X X Y 2 , X X Z 2 , Y X X 2 , 
Y{Y 2 , or Y\Z 2 , which is again 1/2. It is not hard to 
see that this pattern holds for any choice of Hp — the 
probability that it commutes with a randomly chosen Uj 
is always 1/2. Thus, the probability that all m of the Uj 
commute with Hp is simply l/2 m . 

Now impose the constraint that we wish to keep a par- 
ticular term H a while eliminating Hp in this way. Instead 
of choosing Uj completely randomly, we instead generate 
a random product of Paulis, and then check to see if it 
commutes with our desired term, H a . If it docs, which 
happens with probability i, then we add it to the set 
{Uj}, otherwise we repeat the process. A simple case 
analysis now shows that the probability of a given Uj in 
this set commuting with Hp is |. Thus, the probability 
that all m of the Uj chosen in this fashion commute with 
Hp is again l/2 m . 

In general there will be many terms in the expansion 
of H that we wish to eliminate while isolating H a . The 
probability that _ff( m ) contains a term other than H a 
(i.e., our procedure has failed) is certainly no greater than 
the sum of the probabilities that the procedure failed to 
eliminate each term Hp, that is N/2 m where N is the 
number of terms in H that must be eliminated. □ 



VII. SUMMARY AND FUTURE DIRECTIONS 

We have examined the problem of simulating Hamil- 
tonians using a fixed multi-qubit Hamiltonian and local 
unitary operations. We have provided a classification 



scheme for the simulations that are possible given these 
resources. In particular, we have demonstrated that 
there are only two physically distinct classes of entangling 
Hamiltonians up to local unitary operations. One class, 
the class of odd entangling Hamiltonians, when assisted 
by local unitary operations, can simulate all odd Hamilto- 
nians but nothing else. The other class, the class of non- 
odd entangling Hamiltonians, can simulate all Hamilto- 
nians and are thus dynamically universal. We have also 
demonstrated that all odd entangling Hamiltonians can 
be made universal through the use of a simple encoding 
scheme. Furthermore, the Lie algebras generated by the 
odd entangling Hamiltonians together with local unitary 
operations were shown to be isomorphic to either so(2") 
or sp(2"), depending on whether n is even, or odd, re- 
spectively. 

In this paper we have made a distinction between sets 
of resources that are universal for quantum computa- 
tion and those that are dynamically universal. This 
distinction has been necessary because we can not find 
an efficient means to simulate the Hamiltonians used in 
quantum computation with an arbitrary fixed entangling 
Hamiltonian and local unitaries. We have demonstrated 
that when restricted to fc-local Hamiltonians with lower- 
bounded coefficients dynamic universality is equivalent 
to universality for quantum computation. The resolu- 
tion of when this equivalence holds in general would be 
an interesting contribution to the study of quantum dy- 
namics. 
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